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We have performed a systematic angle-resolved photoemission study of chemically substituted 
cuprates Bi2Sr2-c C Bi cc Cu06+<5. We observed that the Fermi surface area shrinks linearly with Bi sub- 
stitution content x, reflecting the electron doping nature of this chemical substitution. In addition, 
the spectral linewidth broadens rapidly with increasing x, and becomes completely incoherent at the 
superconducting-insulating boundary The d-wave-like normal-state gap observed in the lightly un- 
derdoped region gradually evolves into a large soft gap, which suppresses antinodal spectral weight 
linearly in both the excitation energy and temperature. Combining with the bulk resistivity data 
obtained on the same samples, we establish the emergence of the Coulomb gap behavior in the very 
underdoped regime. Our results reveal the dual roles, doping and disorder, of off-plane chemical 
substitutions in high-T c cuprates and elucidate the nature of the quantum electronic states due to 
strong correlation and disorder. 

PACS numbers: 71. 25. He, 74.25. Jb, 74.72. Hs, 79.60.Bm 



Since high temperature superconductivity in cuprates 
is achieved by adding extra carriers into an antiferro- 
magnetic Mott insulator, how the electronic structures, 
such as Fermi surface (FS) and energy gap (including 
the superconducting gap and the normal state pseudo- 
gap), evolve with doping is critical to understanding 
the elusive superconducting mechanism. Direct spectro- 
scopic probe of the electronic structures, such as angle- 
resolved photoemission spectroscopy (ARPES) and scan- 
ning tunneling microscopy (STM), are mostly done on 
the Bi2Sr 2 Ca n Cu n+ i02n+6 compounds, whose doping 
range is limited through oxygenation. More heavily un- 
derdoping can only be achieved by chemical substitution 
[TJ |2]. By substituting trivalent Bi ions for divalent Sr 
ions into the apical plane of Bi 2 Sr2CuO6+5(Bi2201), we 
have further pushed the underdoping threshold contin- 
uously to the insulating phase [3 j. We establish a con- 
crete linear relationship between the reduction of the hole 
concentration and the Bi substitution through ARPES 
by measuring the reduction of the Luttinger volume en- 
closed by the FS. In addition, we observe that quasiparti- 
cle coherence vanishes at the superconducting-insulating 
boundary. More importantly, the d-wave-like pseudo- 
gap observed in the lightly underdoped region gradually 
evolves into a large soft gap, which suppresses antinodal 
spectral weight linearly in both the excitation energy and 
temperature. We propose that the pseudogap in the un- 
derdoped regime evolves into a Coulomb gap beyond the 
superconducting-insulating phase boundary. To support 
this finding, we performed transport measurements on 
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the same samples and found that the resistivity data are 
consistent with hopping transport in the presence of a 
Coulomb gap in the highly underdoped regime, charac- 
teristic of strongly disordered correlated systems. The 
Coulomb energy scale derived from transport agrees re- 
markably well with the gap energy in ARPES spectra. 
The observed universality of this pseudogap evolution re- 
veals the dual roles, underdoping and disorder, of chem- 
ical substitution in high-T c cuprates. 

High quality single-crystals of Bi2Sr2- ;E Bl c Cu06+ ( 5 
(Bi-Bi2201) were grown by the traveling solvent float- 
ing zone method. T c of superconducting samples were 
determined by AC susceptibility, the superconducting 
transition width is about 0.5 - 1.3 K. More technical 
details have been described previously [3 j. ARPES ex- 
periments were performed at the beamline U1NIM of 
the Synchrotron Radiation Center in Wisconsin, and the 
beamline 12.0.1 of the Advanced Light Source in Cali- 
fornia. Energy and momentum resolutions have been set 
to - 10 - 20 meV and - 0.02 A -1 , respectively. All 
samples were cleaved and measured in situ in a vacuum 
better than 8 x 10 -11 Torr on a flat (001) surface. 

We show in Fig. [I] the evolution of the Fermi surface 
of Bi2Sr2- a; Bi :E Cu06+5 at 20 K for six different Bi con- 
tents from x = 0.05 to 0.4. The Fermi vectors extracted 
from the momentum distribution curve (MDC) peaks are 
consistent with the intensity plots and allow a precise 
determination of the underlying FS, which is defined as 
a contour of the minimum gap locus [4]. We report in 
Fig. [l^ all the underlying FS extracted for the vari- 
ous Bi concentrations using an effective tight-binding fit. 
We notice a smooth evolution in the shape of the FS. In 
particular, the FS at the antinodes moves away from the 
M(7r, 0) point as the Bi content increases. More impor- 
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tantly, as shown more clearly in Fig. [T^, the size of the 
Y(7r,7r)-centered holelike FS decreases continuously with 
the Bi concentration, indicating the underdoping nature 
of the Sr 2+ ^Bi 3+ substitution. 
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FIG. 1: FS mapping of Bi 2 Sr 2 -xBi ;r Cu06+(5 at different Bi 
contents, (a)-(f) Plots of ARPES intensity integrated within 
ou = 0±12 meVfor x = 0.05, 0.1, 0.14, 0.2, 0.31, 0.4. The red 
circles are extracted from MDC peaks, and the lines are the 
tight-binding fits to the circles. All the samples are aligned 
along FY parallel to A, as shown by the arrow in (d). (g) 
FS contours from tight-binding fits, (h) Doping effect in- 
duced by the substitution of Bi (dots, the thin line is the 
best linear fit) and La (solid line) measured by Hall measure- 
ments 5 . (i) Corresponding phase diagrams of Bi-Bi2201 and 
La-Bi2201, in comparison to a generic phase diagram of the 
high-T c cuprates [6]. 



We extract experimental values of the effective hole- 
doping p and plot them as a function of the Bi substi- 
tution x in Fig. [TJi, by using Luttinger theorem, which 
states that the volume of an enclosed FS is proportional 
to the carrier concentration 1+p. While the derived dop- 
ing p decreases linearly with x, as shown in Fig. [TJi, the 
linear fit gives p = 0.182 — 0.36x, which indicates that 
each substituting Bi atom removes only 0.36 holelike car- 
rier from the Cu02 plane. A similar result, albeit with 
an even smaller doping efficiency, has been observed in 
Bi2Sr2_ x La ;E Cu06+ ( 5 (La-Bi2201) by Hall measurements 
[5], also plotted in Fig.JTJi as comparison. The reason for 
this surprising behavior is not clear yet. Although one 
would expect that instead of the usual trivalent ions, the 
divalent ions of Bi or La could be present, measurements 
of Bi core levels in our samples are inconsistent with this 
hypothesis. One possibility to explain this phenomenon 
is that the oxygen dopant concentration 5, which is not 
known precisely in these materials, increases with x, com- 
pensating the extra charge of the Bi 3+ ion, and maintain- 



ing total charge neutrality 25 = p + x. Spectroscopic evi- 
dence has led to a similar scenario for the electron-doped 
cuprates, where the Ce 4+ dopant ions tend to form pairs 
with the extraneous oxygen ions [7] . 

In Fig. [TJ, we construct the phase diagram of Bi- 
Bi2201 based on transport /susceptibility measurements 
[3] and the relationship of p vs x displayed in Fig. [TJi. 
This phase diagram, when rescaled along the tempera- 
ture axis, matches well with the one constructed from the 
Hall measurements of the La-Bi2201 [5]. We note that in 
both Bi and La substituted samples, T c vanishes around 
10% hole concentration, which is larger than the 5% crit- 
ical value found in many cuprates [6]. This suggests 
that underdoping is not the only effect introduced by the 
Sr 2+ -^Bi 3+ substitution. Beyond the chemical doping 
necessary to vary the carrier concentration in cuprates, 
a host of experiments have provided evidences that the 
ionic and electronic structures outside the Cu02 planes 
have important effects on the low-energy electronic states 
and the superconducting properties [8j |9j [TUJ [11] . In par- 
ticular, the substitution of Sr 2+ by Ln 3+ (Ln = La, Pr, 
Nd, Sm, Eu, Gd, Bi) in Bi 2 Sr2- cc Ln ;E Cu06+5 leads to 
a critical temperature that significantly depends on ion 
radius mismatch (Ar): at x = 0.4, T c ~ 30K for La- 
substitution which has the smallest Ar, while T c ~ OK 
for Bi-substitution which has the largest Ar [9]. These 
results have been cited as evidence of strong dependence 
of T c to apical site (A-site) disorder. However, our results 
suggest that this is not the complete story. As shown 
in Fig. [TJi, at the fixed substitution level x = 0.4, the 
hole doping levels are different for La-Bi2201 (p ~ 0.16, 
optimally doped) and Bi-Bi2201 (p ~ 0.05, heavily un- 
derdoped). Nevertheless, there is a factor of 3 in terms 
of the maximum T c 's in these two systems, suggesting 
the superconducting properties of Bi2Sr2- ;E Ln ;E Cu06+ ( 5 
are influenced by both charge underdoping and lattice 
disorder, the dual characters associated with chemical 
substitution. 

Another interesting phenomenon we observed is the 
dependence of the quasiparticle (QP) spectral coherence 
on the Bi content. In Fig. [2j we compare the near- nodal 
and antinodal energy distribution curve (EDC) lineshape 
for various Bi contents. Sharp EDC peaks are observed 
in both directions at low Bi substitution levels. The peak 
broadens as x increases and is not observed in the heav- 
ily substituted samples. This loss of coherence is also 
accompanied by an opening of a soft energy gap char- 
acterized by the suppression of spectral weight in the 
vicinity of the Fermi energy (E.p). The crossover be- 
tween the two regimes occurs around x = 0.20, which 
corresponds to a doping of p ~ 0.1, the superconducting- 
nonsuperconducting phase boundary at zero tempera- 
ture. We plot the width of EDCs at kp for the antinodal 
region in the inset of Fig. [2}d. One clearly sees a signifi- 
cant line width broadening when p < 0.1, corresponding 
to x > 0.2. This may suggest that the superconductiv- 
ity is closely correlated to the QP coherence. We also 
note that the opening of this soft gap first appears in 
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FIG. 2: Comparison of EDC lineshape of Bi-Bi2201 at differ- 
ent Bi contents, (a) Near the node, and (b) near the antinode. 
The inset in (a) and (b) with black color show the k posi- 
tions where those EDCs are taken. The inset in red color in 
(b) shows full-width-at-half-maximum (FWHM) of antinodal 
EDCs at Kf. The dashed line is a guide to eyes. 



the antinodal region and spreads out to the nodal one 
upon increasing substitution, as indicated in Fig. [2] A 
similar but smaller nodal gap has been also observed pre- 
viously in several lightly doped high-T c cuprates, and was 
attributed to disorder effect [T2l. 
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FIG. 3: (a) Temperature dependence of symmetrized EDCs 
at Kf in the antinodal region for x = 0.2 sample, and (b) for x 
= 0.4 sample, (c) Zoom- in view of the 20K EDC shown in (b) 
near Ep and a linear fit to it. (d) Temperature dependence 
of EDC intensity at Ep (dots) and the corresponding linear 
fit (line) for x — 0.4 sample. 



Perhaps the most surprising finding from our study 
is the observation of different antinodal pseudogap be- 
haviors at low and high Bi contents. The temperature 
dependence of the symmetrized antinodal EDCs is given 
in Figs. [3^ and b for x = 0.2 and x — 0.4 samples, re- 
spectively. The antinodal gap is filled in and disappears 
above 80 K for the x = 0.2 sample, similar to the con- 
ventional behavior of the pseudogap observed by STM for 
similar samples [13] . The larger gap of the x = 0.4 sam- 



ple behaves differently in the following ways: (1) it lacks 
a well-defined spectral structure to define the value of 
gap, although the suppression of spectral intensity starts 
from ~ 0.2 eV; (2) while its overall spectral shape is lit- 
tle affected within the temperature range up to 280 K, 
the intensity at Ep is linearly proportional to tempera- 
ture, as shown in Fig. |3]i; (3) the spectral intensity is 
suppressed to zero at Ep linearly in energy at low tem- 
perature, as shown in Fig. This reminds us of the 
classical Coulomb gap (CG) behavior in strongly disor- 
dered two-dimensional systems [14 , and suggest that the 
pseudogap in the underdoped regime has evolved into 
a Coulomb gap following the superconductor-insulator 
transition. 

To test this finding, we performed transport measure- 
ments on the same samples to look for evidence of the 
Coulomb gap. The temperature dependence of the in- 
plane resistivity of Bi2Sr2-^Bi iC Cu06+ ( 5 at different Bi 
content x is plotted in Fig. [4] as the natural logarithm 
of resistivity versus T~ 2 . The superconductor- insulator 
transition around x = 0.2 is clearly visible. For x > 0.2, 
the resistivity curves show approximately linear behav- 
ior over a wider range of temperature, consistent with the 
classical hopping resistivity in the presence of a Coulomb 
gap, p(T) = poexp(To/T) 2 , expected for a disordered 
insulating system with long-range Coulomb interaction 
[14]. Here To = e 2 /^£ is the long-range Coulomb en- 
ergy scale determined by the dielectric constant n and 
the localization length £. From the slope of T-dependent 
curves in Fig. [4^i, we obtain T and plot them in Fig. 
(red dots). For the £=0.4 sample, T ~ 400K, which 
explains qualitatively the soft gap in the antinodal spec- 
trum visible even at T = 280K shown in Fig. [3)3. Indeed, 
extrapolating the linear T-dependence of ARPES inten- 
sity at Ep (shown in Fig. [3]i) allows a rough estimate of 
the crossover temperature Tq G ~ 500 K above which the 
suppression of the spectral weight ("gap") is completely 
filled. Thus the crossover temperature {T GG ) obtained 
from ARPES is consistent with To derived from the re- 
sistivity, suggesting that the same Coulomb energy scale 
is involved in both ARPES and transport. 

The ARPES spectra shown in Fig. (3}d allows an esti- 
mate of the magnitude of the energy gap at the antin- 
ode for the £=0.4 sample, A ~ 0.19 eV. This, com- 
bined with the derived value of Tq ~ 400K for the same 
sample (see Fig. [4J3), gives a ratio of A /To ~ 5. It 
is known from the theory of the classical Coulomb gap 
[14], A /To = #o(0 2 ^o where go is unperturbed density 
of states. A reasonable value of go ~ 2.1 state/eV cell 
was provided by band theory for Bi2201 [15 J. Taking 
£ ~ 8 unit cells, we estimate that the ratio predicted 
by the Coulomb gap theory to be ~ 5. This remark- 
able qualitative agreement between the experimental ra- 
tio and theoretical one strongly supports the Coulomb 
gap nature of the observed soft gap in the photoemis- 
sion spectra. It is interesting to note that To becomes 
smaller as the Bi content decreases and becomes zero as 
x becomes 0.2, the superconductor-insulator boundary 
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revealed by the combined effects of underdoping and dis- 
order. We caution here that the conventional Coulomb 
gap usually referred to the linear suppression of the den- 
sity of states (DOS) in an isotropic system, not to the 
spectral function measured near the antinodal region of a 
highly anisotropic material. Nevertheless, the antinodal 
intensity dominates the DOS in the hole-doped cuprates 
due to the van Hove singularity. 
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FIG. 4: (a) Natural logarithm of resistivity versus 
for samples x = 0.1, 0.15, 0.2, 0.31, 0.4 and 
0.5. (b) Lineshape comparison of the antinodal EDCs in 
Bi 2 Sri.6Bio.4Cu06 +( 5, Bis.iSrsCai-.Y^CusOg (from Ref.[R)]) 
and Cai.95Nao.o5Cu02Cl2 (from Ref. 17). (c) Qualitative 
phase diagram of Bi 2 Sr2- ; rBi cc Cu06+(5- 



A natural question to ask is how universal this 
Coulomb gap behavior is in the heavily underdoped 
high-T c cuprates achieved by chemical substitution. 
To answer this question, we compared the antinodal 



EDC of the x = 0.4 sample with heavily underdoped 
Bi2.iSr2Cai_ :c Y :E Cu208 (T c ~ 30K) [16] and nonsuper- 
conducting Cai.95Nao.o5Cu02Cl2 (Na-CCOC) [17] , as 
shown in Fig. (4}d. We found that the lineshape of these 
different samples are very similar, characterized by the 
opening of a large soft gap and the absence of the QP 
peak. 

Interestingly, as shown in the case of 
E^.iS^Cai-^Y^C^Og [16 , the antinodal leading- 
edge gap was found to increase upon more underdoping, 
while the near-nodal gap seems to decrease and be 
proportional to T c . This apparent gap dichotomy (or 
two-gap scenario) has been interpreted as evidence that 
the d-wave pairing gap opens along the FS arc around 
the node while the antinodal gap is a different type 
of gap that may not contribute to superconductivity 
[T6] . However, with the appearance of a Coulomb gap 
in heavily substituted/underdoped samples, we raise 
another possibility, as summarized in the schematic 
phase diagram shown in Fig. [4]:: upon increasing 
underdoping, the d-wave-like leading-edge gap region 
may be influenced or even truncated by the Coulomb 
gap region, which by itself is induced by the intrinsic 
disorder associated with chemical substitution. While 
this Coulomb gap region might push the onset of the 
superconducting phase to a higher doping level, as 
observed in the substituted Bi2201, it may also result 
in a spin glass phase often observed between the anti- 
ferromagnetic insulating phase and the superconducting 
state at low temperatures. 
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